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We study interacting electrons in a periodic potential and a uniform magnetic field B taking 
the spin-orbit interaction into account. We first establish a perturbation expansion for those 
electrons with respect to the Bloch states in zero field. It is shown that the expansion can 
be performed with the zero-field Feynman diagrams of satisfying the momentum and energy 
conservation laws. We thereby clarify the structures of the self-energy and the thermodynamic 
potential in a finite magnetic field. We also provide a prescription of calculating the electronic 
structure in a finite magnetic field within the density functional theory starting from the zero- 
field energy-band structure. On the basis of these formulations, we derive explicit expressions 
for the magnetic susceptibility of B — > at various approximation levels on the interaction, 
particularly within the density functional theory, which include the result of Roth [J. Phys. 
Chem. Solids 23 (1962) 433] as the non-interacting limit. We finally study the de Haas- van 
Alphen oscillation in metals to show that quasiparticles at the Fermi level with the many-body 
effective mass are directly relevant to the phenomenon. The present argument may be more 
transparent than that by Luttinger [Phys. Rev. 121 (1961) 1251] of using the gauge invariance 
and has an advantage that the change of the band structure with the field may be incorporated. 

KEYWORDS: Bloch electrons, magnetic field, correlation effects, magnetic susceptibility, de Haas-van 
Alphen effect, density functional theory 



1. Introduction 

Extensive theoretical studies have been carried out on 
Bloch electrons in a magnetic field B which display many 
exciting phenomena. Among them are investigations on 
fundamental quantities and phenomena such as the ef- 
fective Hamiltonian,^"^ the static magnetic susceptibil- 
j|.y 1,4,5,10-35 g^jj^ ^j^g Haas- van Alphen (dHvA) oscil- 
lation. Since the relevant energy scale is ^bB < 
ImeV with /^b the Bohr magneton, it seems natural to 
start from the energy-band structure in zero field and 
try to include the field effects perturbatively. However, a 
uniform magnetic field in quantum mechanics necessarily 
accompanies a vector potential with a non-periodic lin- 
ear spatial dependence. Hence the procedure is not at all 
an easy task to perform as it apparently looks, particu- 
larly when the electron-electron interaction is taken into 
account. For example, we still do not have a satisfactory 
calculation of the total magnetic susceptibility of metals, 
even for Li and Na. 

The purposes of the present paper are threefold. We 
first establish a definite theoretical prescription to calcu- 
late properties of interacting Bloch electrons in a uniform 
magnetic field on the basis of the energy-band structure 
in zero field. This includes an extension of the density 
functional theory'*^"^'* to a finite magnetic field. We then 
derive explicit expressions of the total magnetic suscep- 
tibility >0) at various approximation levels on the 
electron-electron interaction. Finally, the formulation is 
used to study many-body effects on the dHvA oscillation 
in metals. 

Let us briefly summarize the relevant works on the ef- 
fective Hamiltonian, the density functional theory, the 
magnetic susceptibility and the dHvA oscillation to- 
gether with the extensions considered here. 



We now have a well-established effective Hamiltonian 
in a uniform magnetic field at the one-particle level;^"^'® 
see also Ref. 9 on mathematical aspects. Let {R~R') 
denote the transfer integral in zero field between two unit 
cells specified by R and R' ; it completely determines the 
energy-band structure. Then the transfer integral in a 
finite field is obtained by the replacement: 

t("\R-R') ^c'^'^'^'t(R-R',B), (la) 

where Irr' is the Peierls phase^ and t{R—R', B^O)^ 
t<-°^R-R'). The relation was first obtained by Peierls 
for the tight-binding model without the B dependence in 
t. The extension beyond the tight-binding model is due 
to Luttinger.^ The B dependence in t was taken into 
account by Kohn,^ Blount,'* and Roth.^ As shown by 
Luttinger,^ eq. (la) provides a microscopic justification 
for the procedure: £^°')ik)^£{-iV-^A) with £(o)(fc) 
an energy eigenvalue in zero field, k a wave vector in the 
first Brillouin zone, and A the vector potential, which 
was a key assumption in the Onsager-Lifshitz-Kosevich 
theory for the dHvA oscillation.'^^' '^^ Here, we also con- 
sider interaction effects on the basis of the treatments 
by Luttinger^ and Roth.^ It is thereby shown that the 
self-energy also experiences the change: 

Y.'^°\en,R-R') — >c*^««' S(e„,i2-i?',B), (lb) 

with £„ the Matsubara frequency. We will provide a defi- 
nite prescription of how to calculate S(£„, R—R', B) for 
an arbitrary finite field B. Thus, our approach is more 
powerful than those by the expansion in A^^'^^ that is 
effective only in the zero-field limit. With an application 
to superconductors in mind, the whole formulation will 
be carried out without assuming a specific gauge for the 
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vector potential in such a way that an extension to a 
non-uniform magnetic field may be performed easily. 

The density functional theory (DFT) is regarded now 
as one of the most efficient and reliable methods for the 
quantitative understanding of solids. Its extension to 
a finite magnetic field has been performed by Vignale and 
Rasolt^^ using the current density j{r) as the relevant 
extra variable, and also by Grayce and Harris^^ choosing 
B instead of j{r). In these formulations, however, one 
has to calculate the electronic structure in B from the 
beginning by treating the field and the periodic poten- 
tial on an equal footing. Thus, practical calculations for 
solids have never been carried out. Here, we will provide 
a prescription of obtaining the electronic structure in B 
within DFT based on the known zero-field energy-band 
structure. This two-step procedure may be regarded as 
one of the main advantages of the present formulation 
over the previous ones.^^'^^ We choose the average fiux 
density B = 'V x A as the relevant variable. The total 
moment may be obtained by iW = —{dQ/dB), and the 
external field is found by H = B-iwM/V with V the 
volume. 

The static magnetic susceptibility x(-B^O) has been 
a matter of extensive theoretical investigations,^''''^' 
and we now have several apparently different but es- 
sentially equivalent expressions of x for non-interacting 
Bloch electrons.^' ^^'■'^^'^° Among them, Roth^ gave a 
complete expression including the spin-orbit interaction 
with a clear derivation process. We closely follow her pro- 
cedure to extend the consideration into interacting Bloch 
electrons. It should be noted that % of interacting Bloch 
electrons were already calculated for the orbital part by 
by Fukuyama,^^ Phillipas and McClure,^^ and Fukuyama 
and McClure,^^ and for both the orbital and spin parts 
by Buot^* and Misra et al.^^ However, most of them pro- 
vided only approximate treatments of x- either the ver- 
tex corrections were neglected^''' ^"'"''^ or the frequency 
dependence of the self-energy were not considered.^^'"^^ 
Also, the treatment by Buot^* fails to incorporate vertex 
corrections explicitly in the general expression. We here 
present a complete framework to calculate x various 
approximation levels on the interaction. Particularly, we 
derive an explicit expression of x within DFT. Calcula- 
tions of the susceptibility by DFT have been performed 
only for the spin part,^^"^^'^^'^^'^^ for the orbital part 
neglecting vertex corrections.^^' '^"'■'^ and for the both 
parts but neglecting vertex corrections.^'^ Thus, there is 
no available expression within DFT for the total suscepti- 
bility with vertex corrections. The formula obtained here, 
which incorporates the effects of the spin-orbit interac- 
tion, core polarizations, interband transitions, and ver- 
tex corrections, is expected to enhance our understand- 
ing on the total magnetic susceptibility of solids. It will 
be shown that the vertex corrections of the spin part in 
our formula naturally include the Stoner enhancement 
factor.10'56 

The dHvA oscillation in metals provides unique in- 
formation on the Fermi-surface structures and interac- 
tion effects. The classic theory at the one-particle level is 
due to Lifshitz and Kosevich,^^ who applied Onsager's 
semiclassical quantization scheme^^ to non-interacting 
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band electrons. The interaction effects were consid- 
ered by Luttinger^* and Bychkov and Gor'kov.^® How- 
ever, Bychkov and Gor'kov considered only an isotropic 
Fermi liquid. Also, the work by Luttinger''^ is based 
on a gauge-invariance argument f (°)(fe)-|-E(°)(£„, fe) — > 
£{-iV- ■^A)-|-I](e„, -i'V--^A) without clarifying the 
structure of the self-energy explicitly. Hence it is not clear 
from Luttinger's argument whether it is the momentum 
p or the crystal momentum hk in the first Brillouin zone 
that is really relevant. The replacement procedure for 
the self-energy needs to be established microscopically 
in the same sense as eq. (la) had to. His argument also 
has an ambiguity as to the "non-oscillatory part" of the 
self-energy, which later caused an interpretation^^' that 
the self-energy does not participate in making up the 
quantized energy levels, i.e., one only has to replace the 
one-particle part as f (°)(fe)^£(— iV— ^A). To remove 
the confusion and also to analyze experiments unambigu- 
ously, it will be well worth placing the theory on a firm 
ground. On the basis of the formulation to derive eq. 
(lb), we will present a hopefully clearer argument for 
the many-body effects on the dHvA oscillation. This ar- 
gument also has an advantage that the change of the 
energy-band structure with the field can be taken into 
account. 

This paper is organized as follows. Section 2 provides 
an alternative derivation of the effective one-particle 
Hamiltonian with the spin-orbit interaction. We com- 
bine the advantages in the treatments of Luttinger^ and 
Roth^ to formulate the problem so that extensions (i) 
to interacting systems and (ii) to a non-uniform mag- 
netic field may be performed easily. Section 3 takes the 
electron-electron interaction into account. We establish a 
perturbation expansion for the thermodynamic potential 
and the self-energy in terms of the energy-band structure 
in zero field. It is shown that only the zero-field Feyn- 
man diagrams are necessary. We also construct a DFT 
in a finite magnetic field. Section 4 derives expressions 
of the magnetic susceptibility at various approximation 
levels, including that of DFT. Section 5 studies many- 
body effects on the dHvA oscillation in metals. Section 
6 summarizes the paper. We put A;b = 1 throughout. 

2. Effective Hamiltonian 

2. 1 Hamiltonian and basis functions 

We consider Bloch electrons in a uniform magnetic 
field described by a Hamiltonian with the spin-orbit in- 
teraction and Pauli paramagnetism: 

(2) 

Here m and e(< 0) are the electron mass and charge, 
respectively, P is defined by 

P^P-lA{r), (3) 

with p the momentum operator and A{r) the vector po- 
tential, V(r) denotes the periodic lattice potential, g is 
the electron g factor, s is the dimensionless spin oper- 
ator, B = V X A{r) is the magnetic field, and /x is the 
chemical potential. 
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The eigenfunctions of eq. (2) at A = are the Bloch 
spinors: 



'ipbkir) 



Ubk{r) , 



(4) 



where fe is a wave vector in the first Brillouin zone and b 
denotes a set of quantum numbers for the band and spin. 
They are ortho normal as {tpbk\fpb'k') = Sbb'5kk' and form 
a complete set. For our consideration, however, the wave 
functions (4) may not necessarily be the eigenfunctions 
of eq. (2) a.t A = 0. Thus, we will proceed by assum- 
ing only the completeness and orthonormality of eq. (4). 
Note that, relaxing the conditions in this way, i/'6fc(r) 
may be chosen analytic in k throughout the Brillouin 
zone, as shown by dcs Cloizcaux'''^"''*^ and Ncnciu.'^''''^ 
Those basis functions, which arc analytic in k but not 
necessarily the eigenfunctions of H, were named quasi- 
Bloch functions by des Cloizcaux."''^'"''* 

A set of alternative basis functions was introduced by 
Wannier,^'" which are more suitable for the present pur- 
pose. They are defined as a Fourier transform of V'bfe(^) 
by 

1 V 



WbR{r) 



-ik-R 



i'bkir) , 



(5) 



where R specifies a unit cell and A^c denotes the number 

of unit cells in the system. It hence follows that the Wan- 
nier functions are also complete \ wbR) {wbii \ = 1 and 
orthonormal as (w6h|m'6'h') = hb'^RR' with WbR{r) = 
Wbo{r-R). 

To describe Bloch electrons in a finite magnetic field, 
we adopt the basis functions introduced by Luttinger.^ 
They differ from eq. (5) by only a phase factor as 

^Uir) = e'^-«W6H(r) , 

where IrR is defined by 



(6a) 



'rR 



e 

he 



JR 



dr'. 



(6b) 



with dr' taken along the straight line path from R to 
r. We assume that the functions { wlpf] form a complete 
set, although they are not orthonormal in a finite mag- 
netic field. This latter inconvenience can be removed by 
considering the linear combination: 



(7) 



R' 



and orthonormalize them as {'fibR\'fib'R')=hb'SRR'; this 
will be performed shortly below. 

Wc now transform eq. (2) into a matrix form by using 
eq. (7). Let us introduce the matrices: 



^ = {SbR,b'R') ! 

n' = {{wU\n\<'R'))- 



(8a) 
(8b) 

(8c) 



Note that O is positive-definite Hermitian and reduces 
to a unit matrix as B ^ 0. Now, the orthonormality 
{0bR\0b'R'}=SwSRR' reads 



S}OS = l, 



(9) 



where 1 is the unit matrix: 1_= {Sbb'Snn'). Equation (9) 
is solved easily by choosing S as Hermitian <S^ = 5 as 

5 = 0^1/2 (10) 

We next express the eigcnstatc of eq. (2) as a linear com- 
bination of eq. (7) as tp{r)=J2bR'^b{R')'fibR{'i')- Then the 
eigenvalue problem of eq. (2) is transformed into 



X! '^bR,b'R' Cb'(R') = £Cb{R) , 



VR' 



with 



H = sn's. 



(11) 



(12) 



A couple of comments are in order on eqs. (7) and (10) 
in connection with the treatment by Roth.^ First, our 
basis functions (7) with the phase integral (6b) have an 
advantage over Roth's basis functions that an extension 
to a non-uniform magnetic field can be performed easily. 
Whereas Roth had to assume that the field is uniform 
from the beginning, there is no limitation at this stage 
on the form of A{r) in eq. (6b). Second, eq. (10) is dif- 
ferent from Roth's choice. Indeed, Roth expressed 5 as a 
sum of Hermitian and anti-Hermitian matrices and fixed 
the anti-Hermitian part by imposing that W be semi- 
diagonal up to the order of B"^. However, eq. (10) has the 
advantage that the formulation becomes more transpar- 
ent. It helps to avoid unnecessary complications at the 
one-particle level and make the extension to interacting 
systems easier. Once the orthonormality is endowed as 
eq. (10), the semi-diagonalization may be performed by 
a similarity transformation. It should be noted, however, 
that the physical results are irrelevant of whether one 
carries out such a similarity transformation at an inter- 
mediate step. Indeed, our procedure without further sim- 
ilarity transformations leads exactly the same expression 
for the non-interacting magnetic susceptibility as that of 
Roth. 

2.2 Overlap integral 

We first concentrate on the overlap (Wbj^l'ii'b'R')' where 
the phase factor in eq. (6) yields a line integral from R' 
to R via r. We transform it as follows: 



Irv 



^rR' 



lrr' 



e 

he 



= Irr' + -^{rRXrR' 

he 



I A{r')-dr' 
Jc 

) • / d« / 

^0 Jo 



= IrR' + ■ {tr XTr'). 

Here C denotes closed path R' 
triangle, vr is defined by 



R, 



dvB 

(13) 

R^R! along the 
(14) 



and the second line follows via Stokes' theorem by writing 
the vector area element as dS= (il' — i?)dMX {r — R)<iv. 
The transformation (13) was used previously in deriv- 
ing a transport equation for superconductors with Hall 
terms. Although we have assumed a uniform magnetic 
field here, an extension to a non-uniform field can be 
performed easily as eqs. (31) and (33) of Ref. 61. 
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Following Roth,^ let us introduce the antisymmetric 
tensors: 



(15) 



'2nc " ' ""^ '2m' 
where Ca/j-y (a, 0,^ = x, y, z) is the third-rank completely 
antisymmetric tensor, and summations over the repeated 
index 7 are implied. The quantity ha/s is essentially in- 
verse of the magnetic length squared. Also useful is the 
following identity for an arbitrary function f{rii,rji'): 

f{rR,rR>)wljf^{r)wb'R'{r) 



= 7^ E /(^«' rR')4ki-r)u,'k' (r) e-'=-«+^'='-H' 

" kk' 

= ^Y.^lk{r)u,>k'{r)f{iVk,-iVk')e-''^--+''^'-n' 

kk' 

= -^^[f^-^^k,i^k')ulk{r)uvk'{r) 



xe 



-ik-r R+ik' -r ^1 



(16) 



where we substituted eqs. (4) and (5) and performed par- 
tial integrations over k and k . 

Using eqs. (13)-(16) and the periodicity Ubkir + R) = 
Ubk{r), we obtain an expression for eq. (8b) as 

ObR,b'R' = e'^^^'J exp{ihc,,3r%r^j^,)wlji{r)wb'R'{r)dr 
= ^-^T.^"'-^''~'''^Obb'{k,B), (17) 



where Obb' {k, B) is given by 

Obb-{k,B)= j e''"'^^-<'ulkir)ub'k'ir) 



k'=k 



dr . 

(18a) 

The expression (17) has a general form for any matrix ele- 
ment between eq. (6), i.e., the Peierls phase factor^ q^^rr' 
times a summation over k for the product of Q^l' i^-^ ) 
and some function of k. 

It is useful for practical purposes to expand eq. (18a) 
in terms of B as 



Obb'{k,B) = 5bb' +J2oli}{k,B) 



(18b) 



where o[^^ is the quantity of the order . The terms 
for j = l,2 are given explicitly by 

Oi^ = ih„0{V^Ubk\V^Ub>k) = */i«/3^a;g,„(fe)4,,,(fc) , 

6" 

(19a) 



O^} = -lKpK,p. (Vfe«- Vfe«njfe| Vfe^ V^Ub.k) , (19b) 



7t \V70'V70, 



with 



^bb'ik) = iiubkl'^kUb'k) = -iC^kUbklub'k) 

= J2e"'-''{wbo\r'k\wb'R). (20) 



R 



The second expressions of eqs. (19a) and (20) have been 
derived by noting U6fe(r -|- R)=Ubk{r), reducing the in- 
tegral into a unit cell, and using the completeness and 



orthonormality of {ubk} over the unit cell.^ The last ex- 
pression of eq. (20) is given entirely with respect to the 
Wannier functions of eq. (5). Hence it may be suitable for 
the evaluation of the matrix element x^^, (fc) between the 
core states. Equations (18a) and (19) are exactly those 
obtained by Roth.^ 

A couple of comments are in order before closing the 
section. First, the expansion of eq. (18) is convergent 
if we choose V'&fc(^) as the quasi-Bloch functions of des 
Cloizeaux^'^"^"^^ which are analytic in k. This property 
can be satisfied for a simple band by the eigenfunctions 
of eq. (2), but not for a complex band with band cross- 
ings. This statement holds for every expansion we shall 
encounter below. Second, matrix elements like eqs. (19) 
and (20) cannot be determined uniquely due to the gauge 
degree of freedom: 

W6fc(r-)^e''^'"=U6fc(r). 

However, obsevable quantities such as the magnetic sus- 
ceptibility do not depend on the choice of the phase 4>bk- 
See also the discussions of Roth^ on this point, and those 
by Resta^^ for the electronic polarization. 

2.3 Expression of S_ 

We next derive an expression of S in powers of B start- 
ing from eq. (9) and S} =S. Equation (17) suggests that 
we may expand the matrix element as 

SbR,b'R' = E ^'"-^''-'''^^bb' (fe, B) , (21) 

k 

with SJ {k, B)=S_{k, B). Let us substitute eqs. (17) and 
(21) into eq. (9). We then encounter a summation over 
R" with both the Peierls phase factor e''^««"+*^«"«' and 
the Bloch phase factor e^k (R-R")+ik'-{R"-R') _ rpj^jg g^^^^, 

mation can also be transformed with the procedures of 
eqs. (13) and (16). For example, we obtain 

E ^bR,b''R"^b"R",b'R' 
b"R" 



b"k 



e^HR-R')Sbb„{k)®Ob"b'{k) , (22) 



where the operator 

S(k)(g)0{k) 



is defined by 
^^'^^ S{k')Oik) 



k'=k 



(23) 



Equation (23) is exactly the multiplication theorem of 
Roth.5 

It is worth summarizing the properties of the operator 
(g). First, we realize from eq. (22) that it is associative as 

A{k)(»[B{k)(»C{k)] = [A{k)(»B{k)](8)C{k). (24a) 
Second, it satisfies 

J2Ak)<»B{k) = J2Ak)B{k), (24b) 

k k 

as can be shown with partial integrations over k and the 
antisymmetry ha/3 = —hi3a- 



J. Phys. Soc. Jpn. 



Full Paper 



TaJjafumi Kita and Masao Arai 



By using eqs. (17), (21) and (22), eq. (9) with !S^=!S Note hv" = n°' + iix^H^^^ - n^°^x°'), as can be shown 
is transformed into 

S{k, B) O 0{k, B) ® S{k, B) = l, (25a) « V/3' V^^'Trf \uvk) in eq. (29a) as 
with 1 = ((5b(,'). Let us expand <S in powers of B as 



S{k,B) = l + Y^S^^\k,B). (25b) 

Wc then find from cqs. (18b), (25a), and (25b) that S^^^i 
for j = 1, 2 are given by 

1 



from TTfc = dH^^^ /dk. Next, wc express Ti.'j^^lub'k) and 
I \ub'k) in eq. (29a) as 

b" 

iha'p'Vk ^k \ub'k)= ha'0'^\ub"k) TL^{k)x°'{k) 



b"b' 



5(i) = _^0(i), 

^(2) ^ _i0(2) I £ 0(1)0(1) 

- 2— 8— — ■ 



2.4. Matrix element of Hamiltonian 



We then expand eq. (29a) in powers of B and perform 
straightforward calculations order by order to obtain 
(26a) n'^^^ in eq. (29b). 

The first-order term can be written down easily. Av- 
eraging the resulting expression with its Hermitian con- 
(26b) jugate, we obtain 

2£'(i)= {O(i),H(0)}-/i„^({x",ft£^+7r''}+a"'5) , (31a) 



We next consider Hi 



bR.b'R' 



and w^j^ given by eqs. (2) and (6), respectively. The fol- 
lowing identity is useful for this purpose:^ 



"bfil^Nfc'H') with ^ with 0^^'> given by eq. (19a) and {A,B}=^{AB+BA). 



he nc Jo 

= ^^^")-2^^""«- 



(27) 



It is worth noting that this identity can also be extended 
easily to a non-uniform magnetic field. ^ Using cqs. (13), 
(16) and (27) as well as Ubkir + R) — Ubk{T), the above 
matrix element is transformed into 



Hi 



bR,b'R' 



where W^^, is defined by 
Kb'{k,B) 



>(0) 



Ub'k{r) 



dr . 



(29a) 



k'=k 



with n)^' = n{p^ p+hk,A = Q), TTfc = dnl^'^/dk = 

^ [p+hk-g^^s X VV(r)] , and a°'f^ given in eq. (15). 
We now expand eq. (29a) formally in terms of B as 

00 

Kb, {k,B)= n^} {k) + J2 -X-W {k,B), (29b) 



with 



'H^^,{k)^{ubk\K>\ub,k)- 



(0)U7 



(30a) 



To write down n'l,^{k, B) in eq. (29b) explicitly, we in- 
troduce the following matrices: 



■Kbb'{k) = lubk 



h dk 



eh 



(30b) 



a^b^{k) = {ubk\a''^\ub'k)- 



Ub'k^ 

(30c) 
(30d) 



To obtain H'^^^ compactly, on the other hand, we 
carry out partial integrations using the antisymmetry 
hap = —hfSa- For example, one of the relevant terms are 
transformed as 



,^3ife.(fi-H')7^,^,(fe,B), (28) ^ith 



^ - V^({K'Kubk\ub"k)n4b, 

(,// l L 



+ ^(Vfe"'Vfe««6fe|tl6"fe)VfvM6'(fe) 



{yk'ykUbkiub'k) = iv^'x^b' -Y.^w'4'b' 

b" 

= iVkXbb' — ''^Xbb"Xb'ib' , 



6" 



and /.«,.Vf [(V^".<,;0^^f"J =/^«/3(Vfe"a:^5''0Vf Vf 



b"b' 



Similar calculations are performed for the other terms to 
make H'^^^ compact and symmetric. We finally perform 
the averaging with the Hermitian conjugate. We thereby 

obtain 



H'(2) = {o(2),H(°)} + /,„^/,„,^,<jvf 



+^ix"]T!^' x"' +x°'' n^' x") + 



(31b) 



with {A,B} = ^{AB + BA) and [A,B] = AB-BA. 
Equation (31b) is in agreement with eq. (57) of Roth.^ 

2.5 Effective Hamiltonian 

We are ready to derive an effective one-particle Hamil- 
tonian in a finite magnetic field. Let us substitute eqs. 
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(21) and (28) into eq. (12). Using eq. (22) twice, we ob- 
tain 



R,b'R' — 



Nr. 



with 



(32) 



(33a) 



n{k, B) = S{k, B)(g>]i!{k, B)®S{k, B) 
Let us expand eq. (33a) in powers of B as 

oo 

n{k, B) = n^°^ (fe) + K^'^ {k,B), (33b) 

with ]iJ'^\k) given by eq. (30a). The expressions of TiS-'^ 
can be obtained order by order by substituting eqs. (25b), 
(29b) and (33b) into eq. (33a) and using cqs. (26) and 
(31). The first-order term is obtained from H^^^ =n'^^^+ 
2{5(i),2i^°)} as 

7£(i)(fc) = -h^f}{{x",hv^ + n^}+a"^) , (34a) 

with {A,B} = ^{AB + BA). The second-order contri- 
bution is calculated from W^^^ = W'^^^ + 2{S(^\ } + 
2{5(i),2£'^i^} + 5«W^°)5« -z/i„^[Vf5(i),Vfe«W^°)] as 



{0(i),2iW} + i[0(i),[0W,2£W]] 
J/ia/3Vf[0W,fir], (34b) 



and use the following identity proved in Appendix II of 
Ref. 2: 

exp(i/R,K_H") exp {-R"-Vr) = exp {-R"-dR) , 
with Or defined by 

Equation (36) is thereby transformed into 

J2 'Xbb' {-idR, B) cv [R) = £cb{R) , (38) 

b' 

where Tibb' r, B) is an operator symmetrized with 

respect to Or as 



Ki-idR, B) = Yt{R", B) e'"-"- . 



(39) 



R" 



Equation (38) extends the well-known rule £{k) 
£{—i9R) in a finite magnetic field^'^^ to include the 
change of the energy-band structures in B. 

Next, Dyson's equation corresponding to eq. (38) is 
given by 

'Yyi'^nh'b"-'Hb'b"{-idR',B)\gb"R',bR{£n) = Sb'bSR'R, 

b" 

(40) 

where gb" r' ,hR{^n) is the non-interacting Matsubara 
Green's function with £„ = {2n + l)iTT {n = 0, ±1, • • • ). 
Suggested by eq. (36), we write the Green's function as 



9b'R',bR{£n) = e'-^^'^ 9b'R',bR{Sn) , 



(41a) 



and substitute it into eq. (40). Using cq. (27), we then 
find an equation for g which is also given as eq. (40) with 
a replacement: 

dR' — > = ^ - i-^Bx(R'-R). 
^ ^ dR' 2hc ^ ' 

It hence follows that gwR' ,bR{sn) depends only on R'—R 

and can be expanded as 



with {A,B}^^{AB + BA) and [A,B]^AB-BA. gb'R'M^n) = ^Y.^^'^'^''' '""^ Sb'b{en,k) . (41b) 



2.6 Schrddinger equation and Green's function 

We finally write down an effective Schrodinger equa- 
tion and the corresponding Green's function for non- 
interacting Bloch electrons in a magnetic field. To this 
end, it is useful to introduce the transfer integral in terms 
of eq. (33) as 



t(il,B)^^^e^'=-«K(fe,B). 



(35) 



Substituting eq. (32) and using eq. (35), the Schrodinger 
equation (11) is transformed into 

Y c''^"-^' tbb' {R-R', B) cb' {R') = £cb{R) . (36) 

b'R' 

Thus, the Hamiltonian is given as a product of the trans- 
fer integral and the Peierls phase factor. Note that the 
transfer integral also depends on the magnetic field here. 
We can provide an alternative expression to eq. (36). Let 
us change the summation in the above equation from R! 
to R!' = R-R', express Cb'{R-R") = e--^" '^«C6'(i2), 



Let us substitute eq. (41b) into the equation for g, ex- 
press {R'-Ry''<^'-^^ = -iVke'''<^'-^\ and perform 
partial integrations over k. We thereby obtain an equa- 
tion for (?(£„, k) as 

Y ^kk' [i£nl -H{k,B)] g{en ,k') = l, (42) 
fe' 

where the operator k is defined by 

K = k-i-^BxVk'. (43) 
2hc " ^ ' 

Note that cq. (42) may be written alternatively by using 
the operator of eq. (23) as [ie„l — Ti,{k, B)] (g) Sf(e„, fe) = 
1. This latter result may have been obtained directly by 
writing Dyson's equation (40) in terms of the Hamilto- 
nian in cq. (36), substituting eqs. (35) and (41) into it, 
and using eq. (22). 
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3. Interaction Effects 

We now include the two-body interaction: 

U{r-r') = ^^llqe"'^'-'''^ (44) 

into our consideration, where V is the volume of the sys- 
tem. 

The total Hamiltonian Titot is given in second quanti- 
zation by using the basis function (7) and the notation 
v = bR as 
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(45) 

Here c^, is the fermion operator and T-L^u' is given by eq. 
(32). The quantity 1,2; 1,^1,^ is defined by 

(46) 



where S^v' is given by eq. (21) and ^^31,^.^^^^ denotes 



U' 

f3 I/4;f3f4 



= ydrydr'zi;:t(rX,(r)ZY(r-r')<(r')iI?:,^(r'), (47) 
with w^,(r) defined by eq. (6). 

3.1 Bare vertex 

As shown in Appendix A, eq. (46) can be written al- 
ternatively with respect to the Bloch states {■(/'6fc(f")} in 
zero field as 



^^Rl R'^ '^^^ R2 R^ 



VIU2\V[V'2 



y^^^ik-iRi-R[)+ik'-(R2-R'2. 



kk' 



xAt,b',ik\-q,B). (48) 

Here Irr' and Uq are given in eqs. (6b) and (44), respec- 
tively. The quantity A(fc, q, B) may be called a vertex 
function which is defined in terms of the operator ® in 
eq. (23), S in eq. (25) and A' of eq. (A-3) by 



A(fc, q, B) = S{k + q, B)(E)g{k, q, B)®S{k, B) , (49a) 



with k-\-q the wave vector in the first Brillouin zone cor- 
responding to fe + q in the extended zone scheme. Thus, 
k in A(fc, q, B) belongs to the incoming electron, q is an 
additional wave vector from the interaction, and fc + q 
specifies the outgoing electron. In zero field without the 
Peierls phase factors, the interaction (48) in k space may 
be expressed diagrammatically as Fig. 1. We shall see 
below that this diagram in zero field suffices for the dia- 
grammatic calculations of the properties in a finite mag- 
netic field. 

ft is also possible to expand eq. (49a) in powers of B 

as 

00 

A(A;, q, B) = A(") (fe, q)+Y^ A^^^ {k, q, B) . (49b) 




A(k,q) 




A(k:-q) 



Fig. 1. A diagrammatic representation of the interaction (48) in 
the absence of the Peierls phase factors. The vertex A(fc, q) is 
given by eq. (49), and k+q denotes the wave vector in the 
first Brillouin zone corresponding to fc + q in the extended zone 
scheme. 



The expressions of A^-''^ may be obtained with the same 
procedure as that of deriving eq. (34) . The quantity Af"-* 
is given by 

A^l^lik, q) = (^^^le^^'^+o-l^.'fe) , (50a) 



with Kk+q = k + q k + q a reciprocal lattice vec- 
tor. The first-order term is obtained from ^^\k,q) = 
A'(i) (fe, q) + Sj'^ (fc + q) A(") {k,q) + A^o) {k, q)S^^^ (k) to- 
gether with eqs. (26) and (A-4a) as 

A(l)(fe,q) = -^h^0jf(kT^)V^I^°Hk,q) 



vfA(°)(fc,q)l x"(fc). (50b) 



Note that the argument of is k (k + q) when it ap- 
pears to the right (left) of A^*'''(fc, q); keeping this rule in 
mind and dropping the argument of x" , eq. (50b) may 
be written exactly as the hvf contribution of eq. (34a). 
This rule also applies to higher-order terms. 

The second-order term is calculated similarly. We ob- 
tain an expression analogous to the hv contributions of 
eq. (34b) as 

A(^)(fc,q) 



;KpK,p.V^ {^"(fc + 9),2l"'(fc+^)}Vf' A(°\fe, q) 
+ [Vf A(°)(fc, q)] {x" (fc) , (fc) } 
{ai"(Ff^),x"'(fc+^)}vf Vf A("\fc, q) 
f[vfvfAnfe,9)]{2i"(fc),x"'(fc)} 
E"'(fcT^)V;,«A(°\fc,q) 
-Vfc"A("\fc,q)x"'(fc) 
0(1) {k + q)ti^^ {k,q)+ A(i) (fc, q)O(l) (fc) 
(k + q)&^~^ (k + q)A^"^ (fc, q) 



'-ha/sha'fS' 
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-20(i) (fc+g)A(") (fe, g)o(i) (fc) 
+A(°)(fc,q)0(i'(fc)0(i)(fc)' 

^a/jVf {o(i)(feT^)V^"A(")(fe,q) 



(50c) 



5.^ Perturbation expansion 

With these preliminaries, we proceed to the calculation 
of the thermodynamic potential ft = — TlnTre"^'"'/^, 
the self-energy Sj/jy' and Green's function Gi,'i, in the 
framework of the conserving approximation.^"*' 

Let us define the Matsubara Green's function by 



'Tr.c,,(r)4)e'' 



(51) 



which as Z^^O reduces to g^'ui^n) of eq. (40). As shown 
by Luttinger and Ward,^'^ can be written as a func- 
tional of G as 

f7 = -T^ Tr{ln[2£ + S - ie„l] + S G}e*""°+ + $[G] . 

n 

(52) 

Here 0+ is an infinitesimal positive constant and $ de- 
notes contributions of all the skeleton diagrams in the 
bare perturbation expansion for f2 with G used as the 
propagator. The self-energy is obtained from $ by 

1 (5$ 

5^i/iy'(£n) = 77; — j—T- (53a) 

With this relation, f2 is stationary with respect to a vari- 
ation in G satisfying Dyson's equation: 

[ienl ~n~ S(e„)] G(e„) = 1 . (53b) 

The conserving approximation denotes approximating $ 
by some selected diagrams and determining G and E self- 
consistently by eq. (53).^"*'^^ Its advantages may be sum- 
marized as follows: (i) both the equilibrium and dynami- 
cal properties can be described within the same approxi- 
mation with various conservation laws automatically sat- 
isfied; (ii) the Laudau Fermi- liquid corrections,^^' or 
vertex corrections in a different terminology, are auto- 
matically included. For example, the lowest-order con- 
serving approximation is nothing but the Hartree-Fock 
approximation where and S|/Ji are given by 

^ = Z-/i/ji/2;i/J,y^ [Gi/Ji/j (Snj )G,y^l/2 (e^a) 

-G,;,,(e„JG,^,,(£„J]e-"i''++-"=''+ , (54a) 

(54b) 

respectively. 

It follows from eqs. (32) and (48) that Eiy^/(e„) and 
(e„) can be expanded as 



^g./c.(H-R')s,,,(e„,fc), (55a) 



KY 



< l < 





Fig. 2. Some typical diagrams for the self-energy: (a) the Hartree 
and Fock terms; (b) the second-order exchange contribution. 



5]e''=-(^'-^^Gfc-fc(e„,fc). (55b) 



This is proved by induction as follows: First, the non- 
interacting Green's function gu'vi^n) is already given in 
the form of eq. (55b) as eq. (41). We next assume the 
expression (55b) and substitute it into eq. (53a) to cal- 
culate Ej,^/ (e„) order by order by using eq. (22). We then 
find that the self-energy can also be written as eq. (55a), 
i.e., the same expression as the non- interacting Hamilto- 
nian (32). It hence follows that Gi/v(e„) may be written 
as eq. (55b). This completes the proof. 

To see how to calculate E(£„, k) practically, we present 
expressions of E^^(A;) obtained from eq. (54b), see also 
Fig. 2(a), and E^''(e„,fe) which corresponds to the dia- 
gram of Fig. 2(b): 



V ^ 



Y^l^KMk. K)Y,'^TK{k\ ~K)G{en, k') 



K 



-^UqA{k + q, -g)(g)G(e„, k + q)(»A{k, q) 



(56a) 



E2e(e„,fe) ^I-Y^Y.UqUq,[K{k+q' ,-q') 



V 



nin2 qq' 



k + q') 1^ A{k + q+q', -q) 



®G_{£n2,k + q+q') ® Mk + q,q') 



(^G{£m,k + q)® Mk,q)\ , (56b) 

where ® operates on the k' or k dependences of the ad- 
jacent functions. Note that removing the operator ® in 
the above expressions yields the self-energies at B = Q. 
We hence have a simple rule to calculate the fe-space 
self-energy for B > Q using only the zero-field Feynman 
diagrams as follows: (i) Label each of the connected elec- 
tron lines with the same wave vector, like k or k' in the 
above examples, and insert the operator (g) of eq. (23) be- 
tween every adjacent A and G. (ii) We can remove any 
single operator ® within every closed electron loop, as 
in the first term of eq. (56a), because of the absence of 
the Peierls phase factor in the final calculation of using 
eq. (22); this may also be proved directly in k space by 
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using eq. (24b). The rules (i) and (ii) apply also to the 
calculations of and G(£„, k). 

We now provide fe-space expressions of eqs. (52) and 
(53). Noting the rule (ii) above, we realize that $ in eq. 
(52) is given alternatively as a functional of G(£„, fc) as 

$ = ^ - ^ Tr eV„, fe, B) G(£„, k, 5)6-"°+ , 

n £=1 fe 

(57) 

with the contribution of the fth-order skeleton dia- 
grams to S. We hence have 

1 (5$ 



E6fe,(£„,fe,B,{G}) 



(58a) 



T 5Gb-bien,k,B) 

Next, Dyson's equation of eq. (53b) is transformed into 

[ienl-U{k, B)-S(£„, fe, B)](g)G(e„, fe, B) ==1, as shown 
by using eqs. (32), (55), and (22). It may be written 
alternatively as 

J2 ^kk' [ienl-U{fi. B) - E(£„, K, B)] G{en, fe', B) = 1 , 
fe' 

(58b) 

where k is defined by eq. (43). Hence the functional (52) 
is given in fe space by 



-T^Tr 



Trfc ln{[ K(k, B) + S(£„, k, B)-i£„l ]} 



+ $^S(£„,fe,B)G(£„,fe,B) 



+m, (59) 



with Trfc denoting the trace in fe space. It satisfies 
Sfl 

0, 



(60) 



6Gb'b{£n,k,B) 

as shown by using eq. (58). 

Finally, eqs. (57) and (58a) may be expanded with 
respect to the explicit B dependences originating from 
the vertex A and the operator as 

oo 

4>(B, {G}) = $(0) ({G}) + J2 iB,{G}), (61) 
E(£„,fc,B,{G}) = s(°)(£„,fe,{G}) 

CSO 

+ Y,^^'\^n,k,B,{G}), (62) 

where ^^^^ and E*^-'^ are quantities of the order . Note 
that and E(°) are different from those at B = due 
to the implicit B dependences through G. It follows from 
eq. (57) that i>(-^^ is connected with E^-'^ as 



= ^ E ^ E T^^'^'K^n, k, B)G{en, fe, B) e^^"0+ . 

(63) 



2£- 

ne fe 



3.3 Hartree-Fock approximation 

We here concentrate on the Hartree-Fock approxima- 
tion. Since there is no £„ dependence in the sclf-cnergy 
E^^(fe), eqs. (57)-(60) can be simpHfied further by per- 
forming every summation over n with the Fermi distri- 
bution function f{e) = l/(e^/'^-|- 1).*^^ The potentials O 



and $ are thereby transformed into the functionals of 
the occupation number defined by 

n(fe,B) = r^G(£„,fe,B)e"^"°+ . (64) 

n 

Indeed, eq. (59) is now approximated (dropping the B 
dependences for simplicity) by 



- ^ TV n(fe) S«''(fe) + , (65) 



with 

$HF 



= ^ I EE ^icTr[n(fe)A(fe, K)] Tr[n(fe')A(fe', -K)] 



.kk' K 



-^WqTrA(fe-Fq,-qr) (g)n(fe-l-q) ® A(fe,q) (g)n(fe) \. 

kg ^ 

(66) 

The self-energy E^^(fe) of eq. (56a) is obtained from 

by 



(67a) 



It hence follows that 12^^ is stationary with respect to a 
variation in ri(fe) satisfying 



n(fe) 



(67b) 



g[W(fe)+E"'' (fe)]/T ^ ^ 

Equations (67a) and (67b) constitute a closed set of self- 
consistent equations. 

Finally, and E^^ can be expanded with respect to 
the explicit B dependences as eqs. (61) and (62), respec- 
tively. The expressions of $^^(-5) and E^^^-'-' are given in 
Appendix B to clarify some of their fundamental prop- 
erties. 

3.4 Density functional theory 

Extensive theoretical studies have been carried out for 
more than three decades to describe atoms and solids 
quantitatively based on the DFT.^^' It is thereby estab- 
lished now as one of the most efficient and reliable meth- 
ods for the quantitative understanding of solids. Hence it 
is well worth applying the present method to the density 
functional theory. 

We consider the cases where the exchange-correlation 
energy Exc is given as a functional of the spin density: 



n±(r) = Tr i(l ± T3) G(r, r, £„) e 



(68) 



where 1, r^, and G arc 2x2 matrices corresponding to 
the spin degrees of freedom with T3 denoting the third 
Pauli matrix. Let us expand G(r, r, £„) in eq. (68) with 
respect to the basis functions of eq. (7) and transform 
the resulting expression with the procedures of eqs. (13), 
(16) and (22). We thereby obtain 



n±ir) = ^J2^Tn{k) p^{k,r) , 



(69) 
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where n(fe) is given by eq. (64) and p^{k,r) is defined 
by 

p±(fe,r) =^(fe)(g)p'=^(fc,r)(g)^(fe), (70) 

with 

^ k'=k 
(71) 

The factor N^. is introduced in eq. (71) to make ^^^(A;, r) 
finite in the thermodynamic hmit. Note that p^{k,r) in 
eq. (69) is just a matrix clement with no information on 
the occupation number. It hence follows that the DFT 
as a functional of n±{r) can be written alternatively as 
a functional of n(fe). 

The thermodynamic potential of the DFT is given in 
a form similar to eq. (65) of the Hartree-Fock theory as 



-T TV Trk ln{ 1 + e-[2iW+5°'"(«)l/^} 



Trn(fc)S°FT(fc) + $ 



iDFT 



(72) 



The quantity ^^^^ is defined by 

^^^^ = ^^[n{r)] + E^,[n4r)], (73) 

where $^[n(r)] and Exc[na{T)] are the classic Coulomb 
energy and the exchange-correlation energy, respectively, 
with n{r) = n+{r)+n-{r) andcr = ±. The former is given 
explicitly by 

°^ ''n{r)n{r') 



$"[n(r) 



-drdr' , 



(74) 



which is equivalent to the Hartree term in eq. (66) with 
UK = ^T^e'^/K^', the K = term should be removed due 

to the charge neutrality of the system. 

The self-energy S'^^^(fc) is obtained from by 

The condition 6ff^'^ /6nb'b{k) = for the thermody- 
namic equilibrium then yields 

-C^) = c[W(fc)+S°'=-T(fe)]/T _^J_ ' (^^^) 

which is the Kohn-Sham equation'*'* in disguise. Equa- 
tions (75a) and (75b) should be solved self-consistently 
for a given E^c- 

The equivalence between eqs. (59) and (72) for the ex- 
act functional Sxc can be established by an argument of 
using the coupling-constant integral. Suppose we change 
the interaction of eq. (44) as U ^ MA and express the 
corresponding thermodynamic potential as O^. The two 
expressions (59) and (72) satisfy the same differential 
equation^^'^^ dQ^/dX = {'Htnt)/>' and the same initial 
condition $7'^='^ = 0o) where (Wj^t) the thermodynamic 
average of the interaction in cq. (45) and Qq denotes the 
non-interacting thermodynamic potential. It hence fol- 
lows that eqs. (59) and (72) are equivalent. 

It is also possible to expand <i)'^^"'" and E^^""" with re- 
spect to the explicit B dependences as eqs. (61) and (62), 
respectively. The expressions of $^^"^0) and S^^'''^^) are 
given in Appendix C. 



4. Susceptibility 

We here study the susceptibility Xa'a of B ^ based 
on eq. (59). We first establish a general procedure to 
calculate Xa'a for a given functional We then special- 
ize to the Hartree-Fock approximation and the density- 
functional theory to derive explicit expressions of Xa'a 
within those approximations. Unlike the treatments by 
Buot^* and Misra et al.^^ where they jumped directly 
to the expression Xa'a = — BH^^'dBo^^ our consider- 
ation will proceed in two stages by first calculating the 
magnetization Ma = —dfl/dBa and then performing an- 
other differentiation with respect to Ba' ■ This approach 
has an advantage that vertex corrections can be incor- 
porated explicitly in the formula. 

4-1 Exact expression 

To derive a formally exact expression of Xa'a, kt us 
start with classifying various B dependences of the ther- 
modynamic potential Cl in Eq. (59) into three groups. 
The first category is that through k in eq. (43). The sec- 
ond category is the explicit B dependences in ]i,{k, B) 
and S(e„, fc, S, {G}). Those in E originate from the B 
dependence in the vertex A of eq. (49) and the operator 
® of eq. (23); see eq. (56), for example. The third cate- 
gory is the implicit B dependence through G, which is 
responsible for the Fermi-liquid corrections.^^' ®^ When 
calculating M = —dQ./dB, however, we need not con- 
sider this third category because of eq. (60). 

The dependence through k needs a special treatment 
due to the property (24b). As shown in Appendix D, it 
does not mix with the other categories for B— >0, yield- 
ing the Landau- Peierls diamagnetic susceptibility^'^* as 



6 \mc) 



x[vf VfG-i(£„,fe)][v/v,^G-i(£„,fe)]| . (76) 



Aa'ot 



This expression agrees exactly with eq. (4.15) of Buot.^* 
We next consider the second category. Differentiat- 
ing eq. (59) with respect to the explicit B dependences 
in H{k,B), S(e„, fc, B, {G}), and $(B,{G}), we obtain 
the relevant magnetization as 

,.s„o,^TrG(s„,fc,B)^^M 



ML 



n 

dHB,{G}) 

dB„ 



(77) 



We further differentiate eq. (77) with respect to Ba' and 
put B = 0, where we also need to consider the implicit 
B dependence through G. Using eq. (58a), we obtain 



-T^e'^^^+^TrG 



a$(2) 



dBa,dB„ dBa,dBa 

dn^^^ aE^^A dG 



dB, 



(78) 



where 'H^^\ ^^^\ and S*--'-* (j = 1,2) are quantities of 
the order B^ with respect to the explicit B dependences 
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as given by eqs. (34), (61), and (62), respectively. To 
find an expression for dG/dBa' in eq. (78), we differen- 
tiate G~^G= 1 with respect to B^t. We thereby obtain 
a loading-order equation for dG/dBa' given in terms of 
dG~^ / dBa' as 



dG 

dBa' 



= G 



dBa' dBa' dBc,' 



G. 



with G = G{B = 0) on the right-hand side. The last 
term in the bracket of eq. (79) originates from the im- 
plicit B dependence though G, giving rise to the Fcrmi- 
liquid (or vertex) corrections.^^' To find an equation 
for dS^°^ / dBa' , we introduce a vertex function as 



-(0) 

" 616^,6262 



.'(1>2) = 



T2 ^Gb.b,(l)^G^fe(2)' 

with 1 = (£„j, fci). Noting eq. (58a), we then obtain an 
expression for dT^^y dBa' in terms of dG_/ dBa' as 



^Sgi,,(l) 

dBa' 



T 



EE 

6260 "2fe2 



(0) ^ JGb'^bA'2) 
bib'M^''^' dBa' 



(81) 

Equations (79) and (81) form a closed set of equations 
to determine dG/ dBa' ■ Introducing the matrix I by 



6363 



2'bibll,b^b22=-^ ^Gb'^63(l)G'fc;j&l(l)ri3i;^;&26^(l>2) , 



the coupled equations are solved formally as 

dGb'^bA^) 



(82) 



dBa 



E E [a -s-^] 

6262 "2*82 



622 



G(2)f^C^ + C^^G(2) 



dBa 



dBa' 



,(83) 



6'52 



with 1 = {6b'^b'Jbib2^kik2Snin2)- Substituting eq. (83) into 
eq. (78), one can check that the symmetry x'a'a —x'aa' 
is satisfied as required. 

Equation (78) with eqs. (80), (82), and (83) enables 
us to calculate the relevant susceptibility x'a'a once the 
functional $ is given explicitly. The total susceptibility 
is then obtained by 

Xa'a ~ Xa'a + Xa'a > (^4) 

with x^a given by eq. (76). The factor (1 — T)~^ in 
eq. (83) originates from vertex corrections for both the 
orbital and spin parts, which have been derived naturally 
in our treatment. It includes the Stoner enhancement 
factor as the intra-band contribution of the spin part. 

4-2 Hartree-Fock approximation 

We now specialize to the Hartree-Fock approximation. 
Due to the absence of £„ dependence in the self-energy 
S^^(fe), every summation over n can be performed^^ by 
using the Fermi distribution function: 

Thus, it is possible to simpliiy the expression of Xa'a 
further. 



Let us adopt the representation of diagonalizing the 
Hartree-Fock energy in zero field as 



nw{k,B = 0) 



E«^(fc,B = 0) 



5bb'£,bk ■ 



(86) 



Then the Landau-Peierls susceptibility of eq. (76) is 
transformed into 



LP-HF. 



5/(6fc) 



(87) 

with m*p,f}{bk) = fi^/(Vf' Vf^fc)- Thus, only the states 
near the Fermi energy are relevant to the Landau-Peierls 
susceptibility with the effective mass in place of the bare 
electron mass. However, this simple result no longer holds 
in those approximations where the self-energy has £„ 
dependence. Equation (87) may be regarded as an ex- 
tension of the Philippas-McClure result for the electron 
gas^^ to Bloch electrons. 

Next, eq. (78) can be written with respect to n{k) of 
eq. (64) as 



■^Trn(fe) 



dBa' dBa dBa' dBa 



-ETr 



dH'-^\k) 9S"^(i)(fc) 



dBa 



dBa 



dn(k) 

dBa' 



(88) 



The quantity dn/dBa' is obtained from eq. (79) by 
adopting the representation (86) and carrying out the 
summation over n. We thereby arrive at the expression: 



dnb'bik) /(6'fe)-/(6fe) 



dBa 



^b'k-^bk 



dHj^lik) ^ d^^,l''>{k) 



.HF(1), 



dBa 



dBa' 



+ 



d^7b''\k) 

dBa' 



(89) 



To find an expression of dS^^ (fe) / dBa' , we introduce 
a vertex function in terms of eq. (66) as 

52$HF(0) 



-.HF(O) 
" bib[,b2b: 

N, 



,(fei,A;2) = A^c 



' 6nb'^bi{ki)Snb'^b2{k2) 
= '^T.[^KAZ'S''uK)AZ^{k,,-K)-Uk2-k.+K 



K 



X aJ^^, (fe2 , fei - fe2 - K)I^b2b', {kiM-ki+K) 
Noting eq. (67a), we then obtain 



\(0) 



(90) 



dBa' 



1 V- pHF(O) , dUb'^bAkl) 

TVc ^ ^bW,^'''"^) dBa' 
6165^ fel 



(91) 

Equations (89) and (91) form a closed set of equations 
to determine dnb'b{k)/dBa'- Defining the matrix: 



jHF 
b[biki,b'2b2k2 ' 



1 /(6;feJ-/(6 ifei)T^HF(O) /, , N 

^^2) 



-^c ^b[ki—^biki 



(92) 
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6j 61 fei 
/fe'ifci)-/(6ifci^ 



,(93) 



with 1= i6b'b{<)bbi6kki)- 

Equations (88) and (93) enables us to calculate x'^a- 
The total susceptibility is then obtained by 



LP-HF 



(94) 



with Xa^L""^ given by eq. (87). 



4-3 Density functional theory 

From a practical viewpoint, it is perhaps most useful 
at present to derive an expression of the magnetic sus- 
ceptibility within the density functional theory. We hence 
consider it in most detail. Since the self-energy does not 
depend on the Matsubara frequency, however, we can ex- 
actly follow the procedure of the previous Hartree-Fock 
treatment. 

Let us adopt the representation of diagonalizing the 
Kohn-Sham single-particle energy in zero field as 

n'^kk) + s«(°)(fe) + i:tf\k) = dbb'^bk , (95) 

where Ti-H}, and S^^/^-* are given by eqs. (30a), 

(B-2), and (C-9), respectively, with Ubvik) = Sbv f {^bk) ■ 
The vertex function corresponding to eq. (90) is defined 
in terms of eq. (73) by 

'>>.W^^''^^''^) - ^'^Snb^^bAki)6nb'M'',) 



V 



+ W E [pliblil^^^r)F^l'ir)pi^^l^^^^^ 

with pI^?^ and F^^J, given by eqs. (C-3a) and (C-7b), 
respectively. We also introduce a matrix: 

1 /(6'ifci ) - /(6ifci ) ^DFT(O) 



jDFT 



N, 



Cb'ifei —ibiki 



-r;;;:;:k(fci,fc2). 



with / the Fermi distribution function of eq. (85). 



(97) 



Table I. Quantities necessary to calculate the magnetic suscepti- 
bility x^f"^ of the density functional theory given by cq. (103). 
Those with underlines are matrices with respect to the band-spin 
index b. 



Ubk (r) 


WbRir) 


Cbfe 


<(fc) 


m*^, (bk) 




(4) 


(5) 


(95) 


(106a) 


(106b) 


(15) 




7r"(fe) 


a"{k) 




•^DFT(a) 


pDFT(O) 


(20) 


(30c) 


(30d) 


(98) 


(99) 


(96) 


jDFT 


A(o) 


p±(0) 




4"'(r) 




(97) 


(50a) 


(C-3a) 


(101a) 


(101b) 


(C.7b) 



It is also useful to define a couple of quantities in terms 
of eqs. (19a), (34a), (B-5), and (C-11) as 

0(«)(fe) ^ = ie^p^^x_P{k)x^{k) , (98) 



d 



n^^^ (fc) (fe) +S"'^(i) (fe)' 



47re^ ~( 



E Jp''^°Hk,r)F^'J,{r)ni^\r)dr. (99) 

<T,cr'=± 



(100) 



Here v" is the renormalized velocity: 

'''''' = ^'''hm^' 

and ^°\k,K), p^^°^k,r), and F^U^) are given by 
eqs. (50a), (C-3a) and (C-7b), respectively. The quan- 
tities and flj?^ (r) are defined by 

/^E/(^^'^)K"^(fc'^) 



bk 



-{0(«)(fe),AW(fe,ii:)}J , (lOla) 



bk 



-{0^"^(fc),p±^"^(fe,r)} 



bb 



(lOlb) 



with 



A'(")(A;,ii:) = ieap^—x^{k)i^''\k,K)x'<(k) , (102a) 

p'±(«)(fe,r) = ie„0^^x'3(fc)p±(o)(fe,r)^^(fe) . (102b) 

Note A'(") = 9A'(i) /dB^ and = /dB^ with 

A'*-^-* and p'^'^^^ given by eqs. (A-4a) and (C-4a), respec- 
tively. 

With these preliminaries, the susceptibility X^f^ can 
be calculated easily by eq. (94) with the replacement of 
the superscript HF^DFT. It may be written as a sum 
of three contributions as 

X^fl = X^^l-^^^ + x^r:-^^^ + x[?-°"^ . (103) 

The first term denotes the Landau-Pcicrls diamag- 
netism corresponding to eq. (87), given explicitly by 



,LP-DFT_ 

bk 



l(eh\^df{^bk) 



Xa'a 6V2c; ^ d^bk m%,g{bk)m;,^{bk) 



(104a) 



with ml,Ahk) = h'l{y^%^bk). 
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The second term in eq. (103) comes from the last term of eq. (95) as 
in eq. (88), i.e., the second-order perturbation with re- 
spect to TiS^^ +^^'^''^\ Its intra- and inter-band contri- 
butions yield the Pauli and van Vleck paramagnetism, 
respectively. It is calculated as 

■ E E ^M^"' 

bifc'jfci 6262^2 



- _ i 



Nr. 4^ 



'.bk ■ 



(105) 



k 



PvV-DFT 



x[(l-X"'')-1 



5'j6ifei,6^62fe2 " 



^(fel) 

f{£,b'^k2)-f{ib2k2) 



Then w"(fe) can be calculated by either of the two ex- 
pressions: 



6 



b'„k2 



b2k2 



(104b) 



where 2^^^ and TY^^'^f") arc given by cqs. (97) and 
(99), respectively, and [/(C^fc2)-/(62fc2)]/(?b^fc2 -62^=2) 
should be replaced by 9/(^(,2fc2)/^62fe2 for ^6^fe2=62fe2- 
Thus, the factor (l—I^^^)^^ from vertex corrections is 
explicitly included in our formula. 

The third term in eq. (103) is due to the first two 
terms in cq. (88). It is calculated from eq. (34b), (B-6), 
and (C-lOb) as 

ea'/3'7'ea/37^E* 



(2)-DFT 



bk 



{x''',x^}bbvtv^. 



~f{ibk){xP\x^U 



2h' 

m 



bk 



+{0("'Hfc),7£("\fc)} 



hb 



+ J E /fe/^)(6'fe-6/e)(Ol^;^Olf,) + Oi^}0^^^) 



bb' 



E^ 



V 



E /-i"^' 



r)F^'J,{r)ni^'\r)dr. 



(104c) 



Equation (103) with eq. (104) is one of the main re- 
sults of the paper. It extends the result of Roth^ to take 
both the Coulomb and exchange-correlation effects into 
account with appropriate vertex corrections. Indeed, one 
can check that X^^"^ reduces in the non-interacting limit 
to the expression obtained by Roth.^^ The formula en- 
ables us to perform non-empirical calculations of the 
magnetic susceptibility in solids based on the density- 
functional theory. Table I summarizes the quantities nec- 
essary for the calculation of Xa^J'- 

For practical purposes, it may be worth transforming 
eq. (100) and m^,^(6fe) = h^/iV^V^^bk) in eq. (104a) 
into expressions without differentiations with respect to 
k. To this end, let us define a transfer integral in terms 



_!^e-*«t6fli?„. (106a) 



R 



Also l/m^,^(6fe) is obtained by 

1 1 9'6fe_^_lv 



ml 



e ^ thR RjjRfii . 



R 



(106b) 

Together with eq. (20), one may perform a calculation of 
X^fJ without recourse to numerical differentiations in fc 
space. 

Three comments are in order before closing the section. 
First, our susceptibility is defined in terms of the mean 
flux density B by Mq,' = Xa'aBa^ which is different from 
the conventional definition M^' = Xa'a^a of using the 
external field H. However, Xa'a found easily from the 
thermodynamic relation H = B — AttM/V. Indeed, Xa'a 
in the crystallographic coordinates is given by 

Xaa 



H 



Xa 



(107) 



- 4:TrXaa/V ■ 

However, the difference between Xaa and Xaa negli- 
gible for most of the non-magnetic materials. It is also 
worth noting that choosing B as an independent variable 
is more favorable when extending the theory to supercon- 
ductors. Secondly, an expression of the spin suscep- 
tibility within the density functional theory was already 
derived by Vosko and Perdew^** incorporating vertex cor- 
rections, which has been a basis of the number of theo- 
retical works on the spin susceptibility of metals.^^'^^'^® 
Compared with it, however, the present formula has a 
couple of advantages that (i) the orbital contribution 
can be calculated on an equal footing with the spin part 
and (ii) vertex corrections are incorporated explicitly in 
the formula without any further approximations. Third, 
any calculations of the orbital susceptibility based on a 
model Hamiltonian, such as the Hubbard model, com- 
pletely fail to incorporate the orbital contributions to 
Xa^a *ind Xa'Q- This is because the key quantity x"(fc), 
which originates from the change of the basis functions 
by the magnetic field, is necessarily set equal to zero in 
those calculations. Thus, we need to include the field ef- 
fect on the basis functions for any practical calculations 
of the orbital magnetism. 

5. De Haas- Van Alphen Oscillation 

We finally study many-body effects on the dHvA oscil- 
lation in metals, limiting our consideration to clean sys- 
tems without impurity scatterings. A definite advantage 
here over the previous studies^^"**^ is that the structure 
of the thermodynamic potential is known explicitly as 
eq. (59). 

Following the original work by Luttinger and Ward,^'^ 
let us regard eq. (59) as a functional of S(e„, k) instead of 
G{en, k). It is also stationary with respect to a variation 
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in S satisfying Dyson's equation (58b). We then write S 
as a sum of the monotonic and oscillatory parts as 



(108) 



As shown by Luttinger within the Hartree-Fock approx- 
imation^* and discussed more generally by Bychkov and 
Go'kov,^® the oscillatory part S"***^ for the spherical Fermi 
surface is smaller than the monotonic part S° by the 
order {fkOc/l^Y^'^ {i^c- the cyclotron frequency). We as- 
sume that the statement holds up to the infinite order, 
as expected for three-dimentional Fermi surfaces. Let us 
expand il(S) at S^S*^ as 



771" u <^^bib', (fei)(5Sb26' (^2) 



;S°-,(fci)Sg-,,(fc2) + ---, 



(109) 



where the term linear in S"*"^ vanishes due to Dyson's 
equation (58b) for S = S° and G = G°. The second 
term on the right-hand side of eq. (109) is of the order 
(hwc/n)^- On the other hand, the first term also has an 
oscillatory contribution originating from the operator k 
in the logarithm of eq. (59); it is of the order {huJc/ /J,)^^"^ 
relative to the nomotonic part of O, as seen from the re- 
sult of non-interacting systems.^'' It hence follows that 
we may neglect the second contribution in eq. (109). 

The term relevant to the dHvA oscillation in f2(S°) is 
given by 



0° 



-TE Tr Trfc In [H{k, B) + S° (£„ , k , B) - «e„l] 



/OO 
/(£) {ln[2£(K,B) + S°^(£,K,B)-£l] 
-00 



- In [n{K, B) S°^* (£, B) - £l] } 



d£ 
2^ ' 



(110) 



where TP^{e, k, B) = E"(£„ —ie— 0+, k, B), and we 
have used eq. (85) to transform the summation over n 
into an integration on the real energy axis. Since the 
oscillation is due to the states near the Fermi level, we 
only need to consider the region £ ~ in the integral 
where ImE°^(£) is infinitesimal positive definite and may 
be approximated by 0+1. The double trace in eq. (110) 
is thereby simplified into 



0° 



/OO 
f[e){\n[hb{e,K,B)-^i-e + iQ+\ 
-00 



■ ln[/i6(£, K, B) — fjb—e — i0+] } 



d£ 
2m ' 



(111) 



whore hh(e,k,B) denotes a characteristic value of the 
Hermitian matrix Hik, B)+ReE°'^(£, fc, -B)+/zl. 

Next, we diagonalize /if,(£, k, B) by solving the eigen- 
value problem: 

hh{e, K, B)(px{e, k, B) = gx{e, B)(fix{s, k, B) , (112a) 

where Lp\ is an eigenfunction and g\ is its eigenvalue. 
The subscript A is composed of the band-spin index 6, 
the Landau level N together with an additional quan- 
tum number distinguishing its degeneracy, and the wave 



vector kz parallel to the magnetic field S || 2;. We then 
define the energy Ex = Ex {B) as the solution of the equa- 
tion: 



Ex=gx{Ex-n,B). 



(112b) 



Note that the quantity I - ^ 



for the state Ex = fJ. 

is inverse of the discontinuity in the single-particle occu- 
pation at the Fermi level so that it is positive. It hence 
follows that e — gx{£) + fJ, is a monotonically increasing 
function for £^0 and i?A~A*- Equation (111) is thereby 
transformed into 



0° 



/OO 
f{e)eis + ^i-gx{e))ds, 

= -TEln[l + e-(^^-'')/^] , 

A 



(113) 



where 9 is the step function. Equation (113) includes the 
main oscillatory part of the thermodynamic potential. 

Since Ex^ n for the oscillatory part of eq. (113), Ex 
may be calculated accurately by a two-step semiclassical 
quantization scheme as follows. We first determine the 
energy Ebk = Ebk{B) by solving 



Ebk = hb{Ebk-fJ.,k,B) 



(114) 



We then adopt the Onsager-Lifshitz-Kosevich proce- 
dure:^^' ^'^ 



S{Ex) = 27r(iV+7) 



\e\B 



(115) 



with N an integer, 7 a constant of the order 1, and S{Ex) 
denoting the area in k' space perpendicular to B spec- 
ified by Ei,k> < Ex and k'^ — kz- Let us substitute those 
quantized energy levels into eq. (113) and follow the pro- 
cedure of Lifshitz and Kosevich.'^^ We thereby obtain a 
theory of the dHvA oscillation where the Fermi surface 
by Ebk{B) replaces the non-interacting Fermi surface of 
Lifshitz and Kosevich.^'^ 

There may be an alternative way to calculate Ex semi- 
classically. We first determine the eigenvalue gx{s,B) in 
eq. (112a) as a function of e by 



S{e,gx) = 2n{N+j) 



\e\B 



(116) 



where S{e, gx) denotes the area in k' space perpendicular 
to B II z specified by h}){e, k' , B)<gx{e, B) and k'z=kz- 
We then obtain the quasiparticle energy Ex by solving 
eq. (112b). Since the characteristic energy of k is fuoc, 
however, the Fermi-surface structures determined by this 
latter procedure will not differ substantially from those 
of the first one. It hence follows that we may adopt the 
first procedure which is certainly more convenient. 

Three comments are in order before closing the sec- 
tion. First, the operator k appears naturally in the ar- 
gument of S° as well as in 7i of eq. (110). This fact 
shows unambiguously the necessity of considering the 
self-energy to make up the quantized energy levels. Equa- 
tion (110) thus removes the confusion on the many-body 
effects of the dHvA oscillation mentioned in Introduc- 
tion; it supports Luttinger's original argument, refuting 
the quantization procedure without the self-energy. 
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Second, the present theory can also treat changes of the 
energy band structure with B by incorporating the ex- 
phcit B dependences in H and E. Indeed, 'H{k,B) can 
be calculated by eq. (34), and E(e„,fc, B) may be ob- 
tained by the procedure of §3.2 with dropping the os- 
cillatory contribution of G(£„, k, B). This effect was ne- 
glected by Luttinger^^ but can have a substantial im- 
portance, particularly when approaching the magnetic 
breakdown/^ It should be noted that this magnetic 
breakdown is beyond the description of the semiclassical 
quantization procedure and we have to solve eq. (112) 
exactly by taking the relevant multiple bands into ac- 
count. Third, we have used in eq. (Ill) a characteris- 
tic value hb{e,k,B) and replaced fe by the operator k. 
This procedure is well defined for a simple band where 
hb{s,k,B) is analytic in k.^ For a complex band with 
band crossings, however, hb{e, k, B) may not be analytic 
in k so that the use of hb{e,K,B) may fail to describe 
some important effects such as band mixings due to k. 
In this situation, one may be required to use a repre- 
sentation where H{k, B) + ReS°^(e, fc, B) + jxl is ana- 
lytic in k and directly solve the eigenvalue problem for 
2£(k, B) +ReS°^(e, k, B)+ii1. 

6. Summary 

We have constructed a many-body perturbation the- 
ory for Bloch electrons in a magnetic field on the 
basis of the energy band structure in zero field. We 
have thereby clarified the structures of the thermody- 
namic potential and the self-energy in a finite magnetic 
field, and provided a microscopic foundation for the re- 
placement procedure on the self-energy: E(°)(e„,fc) 
S(£„, — iV— -^A). This perturbation theory is then ap- 
plied to obtain explicit expressions of the magnetic sus- 
ceptibility X at various approximation levels on the in- 
teraction. The result for the density functional theory is 
given by cq. (103) together with eq. (104) and Table I. It 
incorporates vertex corrections as well as interband tran- 
sitions and core polarizations. The expression enables us 
non-empirical calculations of x- Thus, it will be useful 
to improve our quantitative understanding of the mag- 
netic susceptibility in solids. Finally, we have presented 
a many-body tlicxiry on the dHvA oscillation in metals 
to show unambiguously that the Fermi surface structure 
with interaction effects in zero field are indeed relevant 
to the phenomenon. 

The present formulation may be extended easily to 
a non-uniform magnetic field. Hence an application to 
superconductors will be fairly straightforward. It is also 
desired to make up a many-body theory on the transport 
phenomena of Bloch electrons in a magnetic field. 
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Appendix A: Derivation of eq. (48) 

We here transform eq. (46) into eq. (48). To start with, 
let us consider eq. (47) and rewrite w'^^{r)w'^, (r) in the 



integral with the procedure of eqs. (13) and (16) as 



Nr. 



3«b 



^ ^hsk^fP^k'.ir) ■ 



The quantity w'J^{r')wl,^{r') may be expressed similarly. 
Equation (47) is thereby transformed into 



ifea • Jis — ifeg • Jig +zfe4 • Ji4 — ife^ • Ji^ 



X E"<» Kb', (fes, fe^, q)Kv, (fe4, fe;, -q) , (A- 1) 



where A'^^^^,{k, k' , q) is defined by 



A',,,{k,k',q)^l 



{-k+k'+q)-r ^ih^fi'VS'^k' ,7! 

xwfe'fe' (r)dr . 



Hkir) 



Note Al,,(k,k',q) oc(5_fc+fe'+q i<: with K a reciprocal lat- 
tice vector. It hence follows that k in A^^, (fe,fc',q) can 

be written alternatively as fe = fe'-|-q with k'+q denoting 
the wave vector in the first Brillouin zone corresponding 
to k'+q in the extended zone scheme. Using this notation 
and noting e*('='+'J) -'* = e*('"'+'') -^, eq. (A-1) is simplified 
into 



IJ' , , 



{R3-R'3)+ik' -(Ri-R'i] 



x^e'<'-(«3-H4)z^^A/^^,(fe^g^S) 
1 

xK,b'Sk',-q,B), 
where A^j,(fe, q, B) is defined by 
A'(,6,(fe,q,B) 

^ j e^^'=+«-e*''«'''^^'^-'«t^,(r)«6,fe(r)|^,^^ dr 



(A-2) 



dr . 



(A- 3a) 



with Kk+q = k+q—k+q. The third line follows from the 
fact that the discrete reciprocal vector K^+q is indepen- 
dent of the infinitesimal changes in k or k+q. The wave 
vector fe in A'(fc, q, B) belongs to the incoming electron, 
q is an additional wave vector from the interaction, and 
k+q specifies the outgoing electron. 

Equation (A- 3a) can be expanded in powers of B as 

oo 

A^,, (fe, q, B) = a(°) (fe, q) + Y, ^ ik, q, B) , (A.3b) 

with I^'^^i given by eq. (50a). The quantities A'^^? can be 
calculated with the same procedure as that for obtaining 
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eqs. (31a) and (31b). To start with, let us rewrite 

b" 

or equivalcntly, 

(^.S+^|e^'^'=+'- = l]A(°l(fe,q)(ub.fe| , 

6" 

where use has been made of the completeness 1 = 

Eh" l^*6"fc+^)('«b"fe+^l = Efa" \ub"k){ub"k\ over the unit 
cell. We then substitute the above identities into eq. 
(A-3a), expand it in powers of B to obtain ^■'\k,q), 
and average the resulting two different expressions. We 
thereby obtain the first-order term as 

A'(i)(fe,q) = ^^0^^\k+^)f^°\k,q) + A(°)(fe,q)0(i)(fc) 
-/i,./3x"(fc+^)VfA(°)(A;,q) 
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contribution to is calculated from eq. (66) as 

= :^EE"^Tr[n(fc) A«(fc,if)] 



VfA(°)(fe,q)l^"(fe) L (A-4a) 



-ha/3 



Note that the arguments of O^^^ and are k (k+q) 

when they appear to the right (left) of I^^\k,q). Keep- 
ing this rule in mind and dropping the arguments of 
O^^^ and x°', we can write eq. (A-4a) exactly as the O^^^ 
and hv^ contributions of eq. (31a). This rule also applies 
to higher-order terms. The second-order term is trans- 
formed in the same way into 

= lQ}'\-k + ^)I^"\k,q) + iA(")(fc,q)0(2)(fe) 

+\Ki}K'0'V^ {^"(fe+9),^"'(fe+g)}vf AW(fe, q) 



+ [Vf A(o)(fe,q)] {^"(fc),^"'(fe)} 



rhaBha'B' 



{£"(fe + 9), (fc+^)} Vf A(°\fc, q) 
+ [Vfe^vf A(0)(fe, q)] {x-{k),x-' (fe)} 
x«'(fe+^)Vfe«AW(fe,q) 



-Vfe"A(°)(fe,q)x"'(fe) 



(A-4b) 



Let us substitute eqs. (21) and (A-2) into cq. (46) 
and carry out the procedure of eq. (22) repeatedly. We 
thereby obtain eq. (48). 

Appendix B: *HF(i) ^^^j 5]HF(i) 

In this Appendix we derive explicit expressions for 
$HFO) (J ^ 2) and S«p(i). First of aU, the Hartree 



fefe' K 



xTr[n{k')/^'^\k',-K)] 
= -K0j2Trn{k){x'^ik),Vi^''^°\k)} , (B-l) 

where we have used cq. (50b) to obtain the second ex- 
pression with (fc) given by 

E"(o) {k) = ^Yl ^)Tr[ri(fc') A^") (fc', -K)] . 

k'K 

(B.2) 

On the other hand, the Fock contribution has extra terms 
from the operator (g) as 

$F(i) 

= - "-jTr n(fc) {a(i) (fe+^, -9)n(fe+^) A(0) (fc, q) 

kq 

+^i/ia/3[Vfc«AW(fe+^, -g)]Vf [n(fe+^)AW(fe, q)] 
+ ii/^a/j A(°) (fcT^, -q)[Vfe"n(fe+^)] Vf A(o) (fc, q) } 
= -ha^Y."^ [n(fc){^"(fe), Vfe^E^W(fc)} 



+Mn{k)]{xf'{k),E^^'Hk)} 
+ix"(fc)n(fc)^^(fc)S''(°)(fe) 



(B-3) 



The last expression has been obtained through a calcu- 
lation of using the symmetry /i^^ = —h^a and partial 
integrations with 

EF(0)(fc) = _l^z^,A(o)(feT^, -q)n(fe+^)A(°)(fc,q) . 

(B-4) 

The terms with vfs"(°)(fc) and Vfs^(°)(fc) in eqs. 
(B-l) and (B-3) have the effect of renormalizing the non- 
interacting velocity v in eq. (34a) into j; = Vfc[21(fc) + 
E^^(fc)]/?i. However, there also appears extra contribu- 
tions in which are not directly connected with the 
renormalization effect. 

The corresponding self-energy is obtained easily by 
S"P(i)(fe) = 5$™(i)/fe(fc). The Hartree part is calcu- 
lated from eq. (B-l) as 

S«(i)(fc) = -/i„4^«(fc), Vfs«W(fe)} 

+7 T.^K^^°\k, K)Y^ TVn(fe')A(i)(fe', -K) . 



K 



k' 



(B-5) 

Thus, the Hartree self-energy already has a term which 
can not be expressed in terms of VjS^'''^(fe). The Fock 
self-energy may be obtained similarly from eq. (B-3). It 
also has extra terms besides the one with V^S^*^"-* (fc). 

The second-order functional $^^(2) may be calculated 
similarly. The Hartree part is obtained from eqs. (50) 
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and (66) as 

$H(2) 

= ^Trn(fe) L/3V^, (vf ^«'}, Vf S«W} 



2V 

Kkk' 

(B-6) 

The terms with S^^*^^ (fe) in the above expression have the 
effect of rcnormalizing the non-interacting energy Ti'-'^-' 
and the velocity hv_. Indeed, we can find the correspon- 
dents to them in eq. (34b). On the other hand, the last 
term cannot be expressed as the renormalization effect. 
The Fock part may be calculated similarly. 

Appendix C: *dftO) yP^t(i) 

We here derive expressions 

of $DFTO) for j = 1,2 and 
j^DFT(i)_ Among the two contributions in eq. (73), the 
Hartrce part has already been treated in Appendix B. 
The quantities ^^^^^ and $^(2) .^-^^ given by eqs. (B-1) 

and (B-6), respectively, and E^'^' is obtained as eq. (B-5) 
withZYj^ = 47reVif2 (k^q). 

As for the exchange-correlation part, we here consider 
the cases where i^xc is given by^° 

i^xcK(r),Vn<,(r),An,(r)]. (C-la) 

Let us express na{r) as eq. (69) and regard E^^ ^ 
functional of n{k) instead of na{r). We then expand 
E^c with respect to the explicit B dependences through 
p^{k,r,B) as 

oo 

E^, = E(^J+J2Eg^. (Clb) 

To obtain the expressions of E^c , we expand p^{k, r, B) 
of eq. (70) in powers of B as 

oo 

p± (fe, r, B) = p± W (fe, r) + J2 P^^'^ r, B) . (C-2) 

i=i 

The expansion coefficients are found easily as 



p^J°)(fe, r) ^ nlkivY-^il ± T3) uyk{r) , (C-3a) 

p±(i) = p'±(i)+2{5(i\p±(°)}, (C.3b) 

= ^' ±(2) + 2 {5(2) , p± (0) } + 2 {5 W , p' ± W } 
+ - iKp [V^5(^) , Vgp±(°)] , (C.3c) 



Nr. 
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with obtained from eq. (71) as 



p[ty\k,r,B) EE ih^^mul^ir)] ^(1 ±r3)V^^b,fc(r) , 

(C-4a) 

p',t^'\k,r,B) EE -i/ic./3/ia'/3'[V^Vg'4fe(r)] ^(lirg) 

xV^vft?6.fe(r). (C-4b) 

It follows from eq. (C-1) with eqs. (69) and (C-2) that 

(7) 

Exc for j = 1, 2 are given by 



rt 



E^} = / dr 



5^FW(r)nS^-)(r,B) 



o-=± 



+5,2E^i'Hr,S)i^fj(r)n(V(r,B 



(C-5) 



where na\r,B), Fi^\r), and F^^J,{r) are defined by 
nW)(r,B)= i-^TVn(fe)p-W(fe,r,B), (C-6) 



SEx 



i/=0 



<5[£',n,(r)] ' 



(C-7a) 



^3W = 5^(-l)^-+^'"i5.i5.' 



^[£».n^(r)]5[£>.m^.(»-)] ' 
(C.7b) 

with (D,,p,) = (l,l), (V„-l), iVy,-l), (V^,-l), and 
(A, 1) for z^ = 0, • • • ,4, respectively. 

Noting F^^^ {r + R) = F^^^ (r) , we can further simplify 
the first term in the square bracket of eq. (C-5) by using 
the completeness of {■U6"fc(r)} over the unit cell. To be 
more specific, we expand J2a ^^^\''')^i-'^^^3)^b'k{f) 
in E.^^i'^(Op"^'Hfe,^,-B) as 



(C-8) 



where 



Snb'b{k) 



l-Y.jFi^\r)pll?\k,r)dr 



(C-9) 

denotes the exchange-correlation self-energy in zero field. 
With the same procedure as that of deriving eqs. (31a) 
and (34a), eq. (C-5) for j = l is then transformed into 

= -h.0j2^n{k){x^{k),V^^-<°\k)}. 

k 

(C-lOa) 

Equations (B-1) and (C-lOa) have the effect of turning v 

in eq. (34a) into the renormalized velocity v. 

(2) 

The second-order term Ei^c can be transformed simi- 
larly with the procedure of deriving eqs. (31b) and (34b). 
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(2) 

We thereby obtain an expression of -Bic as 



+ ^ |nW(r,B)F(^j(r)n«(r,B)dr, (C-lOb) 

(7(7' 

which is analogous with eq. (B-6). The terms with 
^xc(o) ^^-j above expression have the effect of renor- 

malizing the non-interacting energy TY^"-* and velocity hv_. 
Indeed, we can find the correspondents in eq. (34b). On 
the other hand, the last term cannot be expressed as the 
renormalization effect. 

The first-order self-energy is obtained from S^^/^"* (fc) = 
as 

S^^(i)(fc) = -/i„0{x"(fc), Vfs'=^(°)(fc)} 

+ 7^E /p'^^°'(fe,r)Ff;,(r)n(V(r,B)dr. (C-ll) 

crcr' 

Thus, S'''^^^^(fe) also has a term which can not be ex- 
pressed in terms of VHE^'^^^k). 

Appendix D: Derivation of eq. (76) 

We here derive eq. (76) valid for B^O by expanding 
the logarithmic term in eq. (59) with respect to the B 
dependence in k up to the second order. Our method is a 
slight modification of the one developed by Sondheimer 
and Wilson^^ and refined by Roth.^ It enables us to treat 
the correlation effects exactly, as seen below. 

Let us define the matrix H_{en, k) by 

H{en,k)=U{k)+U^n,k)-ienl. (D-1) 

It is connected with the Green's fimction as H_= —G_~^ 
and appears in the logarithm of eq. (59) as H_{en, *«). Let 
us denote the eigenvalues of H_{en,k) and H_{en,K.) as 
Ei and E\, respectively. Specifically, £ is given hy £ = bk, 
and A is composed of the band-spin index b, the Landau 
level A'' accompanied by an additional quantum number 
distinguishing its degenerate states, and the wave vector 
kz perpendicular to the magnetic field. Both E£ = Ei{sn) 
and Ex = Ex{en) are complex in general and can be writ- 
ten as 



Ex{en) = E'xien) - iE'l{sn) 



(D-2) 



for example, where E'^ and E" arc some real numbers. 
We will proceed by assuming that E'^ has the same sign 
as Sn, i.e., E'^{sn) = \E'^{en)\sgn{en), which is expected 
from the analytic properties of Green's function. 



Let us concentrate on the case £„ > where E'^ > 0. 
Then the trace of the logarithmic term in eq. (59) is 
transformed as 

5 = TVTrfc ln[^(e„,K)] =^ln£A 

= -E 



-il-t 



■dt 



-TrTrfc 



L 



-il t 



t 



■dt, (D-3) 



with 1_ = 1 — «0+. The final expression enables us to 
expand e~'— in powers of B in k. To this end, we 

adopt the procedure of the perturbation expansion in 
the field theory and define U_{k, t) through 

It can be written explicitly as 



(D.4) 



Uik,t) =Texp 



H'{k,t')dt' 



(D.5) 



where T is the time-ordering operator and Il'{k,t) 

k' 



:^'«(fe)+^'(2)(fc)+... . 



(D-6) 



We perform the expansion of eq. (D-5) up to the sec- 
ond order in B. As noted by Roth,^ the term ff'^"'^-' of 
eq. (D-6) does not contribute since it necessarily yields 
hapi^kH){V^H)=0. Thus, the B dependence through 
K yields no terms first-order in B. It hence follows that 
this dependence of the first category can be treated in- 
dependently from the others in obtaining the expression 
for the zero-field susceptibility. Thus, the expansion of 
U.{k, t) up to the second order is given by 

U{k, t) = l-i f e'^('=)*'ff'(2)(fe)e-'^('=)*'dt' 



-i^^'k^'kH) 



(vf^)(vf^) 



(D.7) 



Let us substitute eqs. (D-4) and (D-7) into eq. (D-3). We 
then carry out a similarity transformation in the trace 
to express S in terms of E^. We also use the identity 
(p = 0,l,2,---): 

(-ii)Pe-'^^* -e"'^-* dP 



-/ 

Jo 



dt = 



= (5polnS,-(l-(5po)(p-l)!GP, (D-8) 

where = —EJ^ is Green's function. We finally 
perform an integration by parts to the contribution 
from the second term in the square bracket of eq. 
(D-7) as T.ki^km^^E,)2G\ = Hki^kEeWkG} = 
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" Sfc(^fe ^fe We thereby obtain a simple expres- 

sion for eq. (D-3) as 



(D.9) 

The case e„ < can be treated similarly. Indeed, we 
only need to change the range of integration in eq. (D-3) 
into (—00,0), replace 1_ by 1+, and remove the minus 
signs in front of the integrals. We thereby arrive at the 
same expression as eq. (D-9). Differentiating the second 
contribution of eq. (D-9) twice with respect to B and 
going back to the original k representation, we arrive at 
eq. (76). 
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